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Abstract 

The relationship between point vortex dynamics and the properties of polynomials with roots at the vortex posi- 
tions is discussed. Classical polynomials, such as the Hermite polynomials, have roots that describe the equilibria 
of identical vortices on the line. Stationary and uniformly translating vortex configurations with vortices of the same 
strength but positive or negative orientation are given by the zeros of the Adler-Moser polynomials, which arise in 
the description of rational solutions of the Korteweg-de Vries equation. For quadrupole background flow, vortex 
configurations are given by the zeros of polynomials expressed as wronskians of Hermite polynomials. Further new 
solutions are found in this case using the special polynomials arising the in the description of rational solutions of the 
fourth Painleve equation. 

1 Introduction 

Vortices are some of the most visually appealing phenomena of fluid mechanics. Examples include the bathtub vortex, 
vortices resulting from flow past obstacles, wingtip vortices, as well as hurricanes and tornadoes. Vortical flows play 
an essential role in many applications and hence it is no surprise that vortex dynamics is an important sub-topic in 
fluid mechanics. A polynomial is a basic mathematical structure, perhaps the simplest and most fundamental. There 
are many situations where approximations of mathematical functions are made using polynomials, for example in 
numerical methods. In this paper we are concerned with the relationship between vortex dynamics and the roots of 
some special polynomials. 

The point vortex equations may be thought of as a discretization of the equations for two-dimensional flow, the 
Euler equations, which is useful for analytical and numerical approximations. Considering the flow plane of the 
two-dimensional fluid to be the complex plane, i.e. the point (x, y) is viewed as the complex number z = x + iy, 
one associates a set of point vortices with a polynomial that has roots at the locations of the vortices and uses this 
representation to study the vortex configurations. 

In two recent papers, Aref [8, 9] reviews the connection between point vortex dynamics and polynomials with 
roots at the positions of the vortices. For stationary vortex configurations the following results have been established: 

1 . N identical vortices on a line are in equilibrium if and only if they are situated at the roots of the classical iVth 
Hermite polynomial, which are related to the Stieltjes relations [50, 51], see equation (3.2); 

2. N identical vortices on a circle are in equilibrium if and only if they are situated at the vertices of a regular 
TV-polygon [37]; 

3. \n(n + 1) positive and |n(n — 1) negative vortices are in equilibrium if and only if they are situated at the roots 
of the Adler-Moser polynomials [10, 12, 32] (see also §2.2 and §3.2.1); 

4. TV positive and N negative vortices form a uniformly translating equilibrium if and only if the vortices of one 
sign are situated at the roots of an Adler-Moser polynomial and the vortices of opposite sign are then at the roots 
of a polynomial of the same degree derived from the Adler-Moser polynomial [10, 12, 32] (see also §3.2.2). 
Further, such equilibria are only possible when N is a triangular number. 

Aref [8, 9] remarks that the relationship between the theory of stationary equilibria and rational solutions of the 
Korteweg-de Vries (KdV) equation 

u t + 6uu x + u xxx = 0, (1.1) 
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one of the original contexts in which the Adler-Moser polynomials arose, is "quite unexpected and very beautiful". 

There has been considerable interest in completely integrable partial differential equations solvable by inverse 
scattering, the soliton equations, since the discovery in 1967 by Gardner, Greene, Kruskal, and Miura [30] of the 
method for solving the initial value problem for the KdV equation (1.1). During the past thirty years or so there 
have been several studies of rational solutions for the soliton equations, applications of which include the description 
of explode-decay waves [40] and vortex solutions of the complex sine-Gordon equation [11, 46]. Airault, McKean, 
and Moser [7] studied the motion of the poles of rational solutions of the KdV equation (1.1) and related the motion 
to an integrable many-body problem, the Calogero-Moser system with constraints. Subsequently Adler and Moser 
[3] expressed the rational solutions of the KdV equation in terms of the Adler-Moser polynomials (see §2.2 below). 
Ablowitz and Satsuma [1] derived rational solutions of the KdV equation (1.1) by finding a long-wave limit of the 
known A^-soliton solutions of these equations. 

The Painleve equations (Pi-Pvi) are six nonlinear ordinary differential equations discovered by Painleve, Gambier 
and their colleagues around the beginning of the 20th century. Their solutions define new transcendental functions as 
they are not expressible in terms of previously known functions such as elementary and elliptic functions or in terms 
of solutions of linear ordinary differential equations and can be thought of as nonlinear analogues of the classical 
special functions [19, 28, 31, 52]. It is well known that Pn-Pyi have rational solutions, algebraic solutions, and 
solutions expressed in terms of the classical special functions (cf. [19, 31]). Ablowitz and Segur [2] demonstrated a 
close relationship between the soliton and Painleve equations. For example, the second Painleve equation (P n ) 

d 2 w , 

— j = 2w A + zw + a, (1.2) 

dz 

where a is an arbitrary constant, arises as a scaling reduction of the KdV equation (1.1) and the fourth Painleve 
equation (Piv) 

d 2 w 1 / dw\ 2 3 A o . 9 , , _ „ N 

10—5- = - — + -w 4 + Azw 3 + 2 z 2 - a )w 2 + 0, (1.3) 
dz 2 2 V dz J 2 

where a and (3 are arbitrary constants, arises as scaling reductions of the Boussinesq equation 

u u + (u 2 )xx ± u xxxx = 0, (1.4) 
and the nonlinear Schrodinger (NLS) equation 

iq t +q xx ±2\q\ 2 q = 0. (1.5) 

Vorob'ev and Yablonskii [56, 57] expressed the rational solutions of Pn (1.2) in terms of polynomials, now known 
as the Yablonskii-Vorob'ev polynomials (see §2.1 below), which are special cases of the Adler-Moser polynomials. 
Okamoto [45] derived analogous polynomials, now known as the Okamoto polynomials, related to some of the ratio- 
nal solutions of Piv (1.3). Okamoto's results were generalized by Noumi and Yamada [43] who expressed rational 
solutions of Piv (1.3) in terms of two sets of polynomials, the generalized Hermite polynomials and the generalized 
Okamoto polynomials (see §2.3 below). 

This paper is organized as follows. In §2 we describe the special polynomials associated with rational solutions of 
the Painleve and soliton equations. In §3 we discuss several examples of equilibrium vortex configurations for vortices 
of the same strength but positive or negative orientation. In §3.2 stationary vortex configurations and uniformly 
translating vortex configurations are shown to be given by the zeros of the Adler-Moser polynomials and related to 
eigenstates of the Schrodinger equation for rational potentials decaying at infinity, which relates to work of Duistermaat 
and Griinbaum [25]. In §3.3, for quadrupole background flow, vortex configurations are shown to be given by the zeros 
of polynomials expressed as wronskians of Hermite polynomials and related to eigenstates of the Schrodinger equation 
for rational potentials with quadratic growth at infinity which relates to work of Oblomkov [44], Veselov [54], and 
Loutsenko [38]. Further using results on the generalized Hermite and generalized Okamoto polynomials associated 
with rational solutions of Piv ( 1 .3) it is shown that there are more solutions in the case of quadrupole background flow 
than might be anticipated from previous work. Finally in §4 we discuss our results. 

2 Special polynomials associated with the Painleve and soliton equations 

2.1 The Yablonskii-Vorob'ev polynomials 

Rational solutions of Pn (1.2), with a = n £ Z, can be expressed in terms of special polynomials now known as 
the Yablonskii-Vorob'ev polynomials, which were introduced by Vorob'ev [56] and Yablonskii [57] and are defined 
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in Definition 2.1. Kajiwara and Ohta [34] expressed rational solutions of Pn (1.2) in terms of Schur functions by 
expressing them in the form of determinants, which is how they are defined below (see also [22, 33, 53]); we remark 
that Flaschka and Newell [27], following the earlier work of Airault [6], expressed the rational solutions of Pn (1.2) 
in terms of determinants. 

Definition 2.1. Consider the polynomials tp n (z) defined by 

oo 

^^„(z)A"-cxp(zA-|A 3 ), (2.1) 
and then the Yablonskii-Vorob'ev polynomial Q n {z) is defined by 

n 

Q n (z) = c n W(<pi,<p 3 ,...,<P2 n -i), c n = Y[(2j + l) n -*, (2.2) 

where W(ipi, <p3, . . . , <p2n-i) is the usual Wronskian. 

Examples of Yablonskii-Vorob'ev polynomials and plots of the locations of their roots in the complex plane 
are given by Clarkson and Mansfield [22], who showed that the roots have a very regular, approximately triangular 
structure. The term "approximate" is used since the patterns are not exact triangles as the roots lie on arcs rather than 
straight lines. 



■ dz, 



2.2 The Adler-Moser polynomials 

Burchnall and Chaundy [3] determined what condition must be satisfied by two polynomials P(z) and Q(z) in such 
that the indefinite integrals 

f P 2 (z) , f Q\z) 

J Q 2 (z) z ' J P 2 (z) 

be rational, i.e. expressible without logarithms. Under the assumption that P(z) and Q(z) have no common or repeated 
factors, the problem is reduced to finding polynomial solutions of the bilinear equation 

d 2 P dPdQ d 2 P n 
Q— T -2—-^ + P— T = 0. (2.3) 

dz dz dz dz 

Airault, McKean, and Moser [7] studied the motion of the poles of rational solutions of the KdV equation (1.1) 
and related the motion to an integrable many-body problem, the Calogero-Moser system, with the Hamiltonian 

N N N 2 

W = E*?+E E t^-^p' 

which are described by the following algebraic system 

N 

E' ( Z .- Zl )3 =°> J = 1.2,--,JV, (2-4) 
fe=i [Z: > Zk > 

N i N 

where ^ means ^ . A remarkable fact, first discovered by Airault, McKean, and Moser [7], is that the system 

(2.4) has no solutions unless = \n{n +1), with n £ Z, i.e. a triangular number, in which case the solutions depend 

on n arbitrary complex parameters. Adler and Moser [3] constructed the polynomials @ n (z) — J\k=i ( z ~ Zfc )' which 
are defined in Definition 2.2 below, whose roots satisfy (2.4) and so all solutions of (2.4) can be obtained. These are 

d 2 

constructed by Adler and Moser [3] using of the Darboux transformations of the operator — , > following the work of 

dz 

Crum [23], though many of the results appear in the paper by Burchnall and Chaundy [13]; see also [39]. 



Definition 2.2. Consider the polynomials 9 n (z) defined by 

oo 

'\ — 

2j 



]T 6 n (z)\ n = exp I zX - f; I , (2.5) 



n=Q \ j=2 
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with Kj 9 j = 2, 3, . . . , arbitrary parameters, and then the Adler-Moser polynomial Q n (z; K n ) is defined by 

n 

Q n {z;K n ) = c n W(6 1 ,e 3 ,...,02 n -i), c n = l[(2j + l) n -*, (2.6) 

where re„ — (k 2 , k 3 , . . . , n n ). The modified Adler-Moser polynomial Q n (z; fi, K n ) is defined by 

G n (z;ti,K n ) = c n e-» z W(9 1 ,6 3 ,...,6 2n - U e» z ), (2.7) 

with fi an arbitrary constant. 

We remark that it can be shown that Adler-Moser and modified Adler-Moser polynomials are related by 

Q n (z;n, k„) = 9„ (z - 1/h;k„) , K n = (k 2 + 1/a* 3 , k 3 + l/^ 5 , ...,K n + l/u 2 ^ 1 ) . 

Further, we note that the Yablonskii-Vorob'ev polynomials Q n {z) are special cases of the Adler-Moser polynomials 
& n (z; K n ) with n 2 = 4 and K n = for n > 3. 

As shown by Burchnall and Chaundy [13], and later by Adler and Moser [3], the Adler-Moser polynomials satisfy 
the bilinear equations 

£|^e.-2^^ + e„ +1 £% = o, (2. 8 , 

dz dz dz dz 

d8 n+ i d6 n _i 2 , 1Q1 

— — — 6 n _i-e„ + i— — — = (2n+l)e n . (2.9) 

In fact equation (2.9) is often used to generate the Adler-Moser polynomials, given <do(z) = 1 and <di(z) = z, with 
the parameters Kj,j = 2, 3, . . . , arising as constants of integration. 

Bartman [12] showed that roots of the consecutive polynomials Q n _i(z) and @ n (z) are the equilibrium coordi- 
nates of \n(n — 1) positive and ^n(n + 1) negative Coulomb charges, with values ±1 respectively, on the plane 
interacting through a two-dimensional (logarithmic Coulomb) potential. Namely, the function 

n(n-l)/2n(n-l)/2 n(n+l)/2 n(n+l)/2 n(n-l)/2n(n+l)/2 

e= y, E ini^--^i+ E E inio-ai+ E E ln i^-Cfei, 

j=l k=j+l 3 = 1 k=j+l j=l fe=l 

has a critical point when Zj, j — 1,2,..., \n{n — 1), and Cfe, k = 1, 2, . . . , |n(n + 1), are roots of 0„_i(z) and 
9„(z) respectively. 



2.3 The generalized Hermite and generalized Okamoto polynomials 

In a comprehensive study of properties of solutions of Piv (1.3), Okamoto [45] introduced two sets of polynomials 
associated with rational solutions of Piv (1.3), analogous to the Yablonskii-Vorob'ev polynomials discussed in §2.1. 
Noumi and Yamada [43] generalized Okamoto's results and introduced the generalized Hermite polynomials, which 
are defined in Definition 2.3, and the generalized Okamoto polynomials, which are defined in Definition 2.4; see also 
[17, 18]. Kajiwara and Ohta [35] expressed the generalized Hermite and generalized Okamoto polynomials in terms 
of Schur functions in the form of determinants, which is how they are defined below. 

Definition 2.3. The generalized Hermite polynomial H m ^ n (z) is defined by 

H m<n (z) = W(H m (z),H m+1 (z),...,H m+n _ 1 (z)) = W({i/ m+j (z)}"r 1 ) , m,n> 1, (2.10) 

with H m fi(z) = H , n (z) = 1 and Hk(z) the kth Hermite polynomial, and has degree deg(H m . n (z)) = ran. 

Examples of generalized Hermite polynomials and plots of the locations of their roots in the complex plane are 
given by Clarkson [17]; see also [18, 20, 21]. The roots take the form ofmxn "rectangles", which are only approxi- 
mate rectangles since the roots lie on arcs rather than straight lines. 

From Lemma 5.7 in [43] it follows that the generalized Hermite polynomials H m , n (z) satisfy the bilinear equations 

[D 2 Z - 2zV z - 2n] H m . n .H m+hn = 0, (2.11a) 
[D 2 Z - 2zV z + 2m] ff ro ,„+i . if m ,„ = 0, (2.11b) 
\pt - 2zD z + 2(m - n)] H m>n+1 .H m+ltn = 0, (2.11c) 
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with D z the Hirota operator defined by 

(2.12) 



jy z (F.G)(z)= ^-{F(z + OG(z-0} 



C=o 

The generalized Hermite polynomial H m , n (z) can be expressed as the multiple integral 

jj.m/2 n™ ^| roo poo n n n 

H m ,n(z) = 9m(m +2»-l 1 )/2 / n WW (x t - xrf ~[[(z - X k ) m CXp (~X 2 k ) d Xl dx 2 . . . dx n , 

J-oo J-oo i=lj=i+1 fe=1 

which arises in random matrix theory [15, 29, 36]. The generalized Hermite polynomials also arise in the theory of 
orthogonal polynomials [16]. 

Definition 2.4. The generalized Okamoto polynomial fl m ,n(z) is defined by 

n ro ,n(*) = W(iTl(z), fl4(*), ■ • • , #3m- 2 (z), # 2 (*), # 5 (z), • • • , # 3n -l(*)), (2.13a) 

= W ({^ 3j+ iW}7 =0 1 , {^ + 2(^)}^: 1 ) , m,n > 1, 
O ro , (z) = W( J ff 1 (z),H 4 (z), . . .,ff 3ro _ 2 («)) = W ({^3 J +i(^)}7=~i 1 ) - ™ > 1. (2.13b) 

= W(H 2 («),ff5(«),..., J ff3n-l(2)) - W ({H 3k+ 2{Z)}"~1) , H > 1, (2.13c) 

with fi ,o(^) = 1 an d H k {z) the fcth Hermite polynomial, and has degree deg(ri m ^„(z)) = m 2 + ra 2 — mw + n. 

We note that the generalized Okamoto polynomial £l m n (z) defined here have been reindexed in comparison 
to the generalized Okamoto polynomial Q m ,n{z) defined in [17, 18] by setting Q m ,n(z) = Cl m + n -i,n-i{z) and 
Q-m,-n(z) = Cl n -i,m+n(z), for to, n > 1. The polynomials introduced by Okamoto [45] are given by Q m {z) = 
^m-i,o(z) and R m {z) = Cl m ,i(z). Further the generalized Okamoto polynomial introduced by Noumi and Yamada 
[43] is given by Q m ^ n (z) ^m— 1 

Examples of generalized Okamoto polynomials and plots of the locations of their roots in the complex plane are 
given by Clarkson [17, 18]. The roots of the polynomial Q m ^ n (z) = £l m + n -i.n-i( z ) with m, n > 1 take the form of 
an m x n "rectangle" with an "equilateral triangle", which have either |to(to — 1) or ^n(n — 1) roots, on each of its 
sides. The roots of the polynomial Q_ m _ n (z) = £l n -i,m+n(z) with to, n > 1 take the form of an to x n "rectangle" 
with an "equilateral triangle", which now have either ^m(m + 1) or \n{n + 1) roots, on each of its sides. Again these 
are only approximate rectangles and equilateral triangles since the roots lie on arcs rather than straight lines. 

From Lemma 5.7 in [43] (see also [38]), it can be shown that the generalized Okamoto polynomials £l m ,n{z) 
satisfy the bilinear equations 

[D 2 - 2zD z + 2(to - 2n + 2)] O m+ i,„ . O m ,„ = 0, (2. 14a) 

[D 2 - 2zD z - 2(2to - n - 1)] Q m ,n+i • On,n - 0, (2.14b) 

[D 2 -2zD z -2(TO + n + 2)] O m>n . Sl m +i,n+i = 0, (2.14c) 

with D z the Hirota operator (2.12). 

The symmetric fourth Painleve (sP IV ) hierarchy introduced by Noumi and Yamada [42], which has the affine Weyl 
group symmetry , for n G N, is given by 



dz 

subject to the constraints 



\ yj+2r-i - ^ yj+2r + 2^- = °> j = 0, 1, . . . , 2n, (2.15a) 



\r=l 



2n 2n 



]T<^- = -2z, ^2fij = h (2.15b) 

3=0 j=0 

where fij, for j = 0, 1, . . . , In, are complex constants and the indices are elements of Z/ (2n+l)Z, i.e. f k +2n+i = <fik- 
We remark that setting n = 1 in (2.15) yields the sPiv system 

^p- + (p (tpi - ip 2 ) + 2fi Q = 0, (2.16a) 
dz 

+ M<P2 - Vo) + 2/ii=0, (2.16b) 

-f- + M<Po-Vi) + 2^2 = 0, (2.16c) 
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with constraints 

<fio + <Pi + ¥2 = -2z, /j, + /Ui + fi2 = 1- (2.16d) 

Eliminating y>i and </? 2 in (2.16) and setting tp = w yields Piv (1.3) with a = fi± — \i 2 and /3 = — 2^q. The sPiv 
system (2.16) was known to Bureau [14], it was subsequently discovered by Adler [4] in the context of periodic chains 
of Backlund transformations, see also Veselov and Shabat [55]. The sPiv system (2.16) has important applications in 
random matrix theory, which are discussed by Forrester and Witte [29]. 

Filipuk and Clarkson [26] studied rational solutions of the sPiv hierarchy (2.15) and expressed them in terms of the 
symmetric Hermite polynomials and symmetric Okamoto polynomials, which respectively are defined in Definitions 
2.5 and 2.6 below and are extensions of the generalized Hermite and Okamoto polynomials. 

Definition 2.5. The symmetric Hermite polynomial Hk 2n (z) is defined by 

H k2n (z) = W ({H kl+J {z)} k ^ , {#fci+fc 2 +fc 3 +.j( z )}jt \ { H k 1 +k2+...+k2 n -i+ 3 ( z )} k j=o *) > ( 2 - 17 ) 
where k 2n = (fci, k 2 , • • • , k 2n ) and Hk(z) is the fcth Hermite polynomial, and has degree 

deg(i?fc 2 „(z)) = k 2n k 2 n-l + (Mn + k 2 n-2)k 2 n-3 + • • • + (k 2n + fc 2 „_ 2 + • • • + k 2 )k 1 . 

Definition 2.6. The symmetric Okamoto polynomial f2fc„ (z) is defined by 

= W {{H (n+ i )j+ i{z)}%\ {H (n+1)j+2 {z)} k ;S \ {H {n+1)j+n (z)}':i~ 1 ) , (2.18) 
where k n — (fei, fe, ■ ■ ■ , k n ) and Hk{z) is the feth Hermite polynomial, and has degree 

n n j—1 n 

deg(ft fe „ (*)) = §n ^ kj {kj + 1) - Yl hk > + Z^ 1 _ ^ ' 

j=l j=2 i=l j=2 

We note that i/fc 2 (z) = Hk lt k 2 ( z ), the generalized Hermite polynomial, and f2fc 2 (z) = 0/^ ,k 2 ( z ), the generalized 
Okamoto polynomial. Filipuk and Clarkson [26] show that pairs of symmetric Hermite polynomials and pairs of 
symmetric Okamoto polynomials satisfy bilinear equations similar to the forms (2.1 1) and (2.14), respectively. 



3 Point Vortex Dynamics 

The equations of motion for N point vortices with circulations Tj at positions zj, j = 1,2,...,N, with no external 
flow field, are 

dz* i „/ r fc 

-r- = 7T^Y — . 3 = 1,2,..., N. (3.1) 

at 2ni f-^ 1 Zj — Zk 

k=l J 

3.1 Identical vortices on a line 

If a vortex configuration rotates as a rigid body with angular velocity fi then 



At 

and so from (3.1) 



dz* 

"' - -in*;, j = 1,2, . . . , n, 



N 

A ^=E'— j = l,2,...,N, 

. ~ Zj z^ 
k=l J 



where A = 27rQ. Suppose that Zj are real, so Zj = z* = Xj, and all the F, are equal, so set F, = 1 for j = 1,2, . . . ,N, 
without loss of generality (by rescaling xj, if necessary), and so we obtain 

N 

Xj = V' , j = l,2,..., N, (3.2) 

which are known as Stieltjes relations [50, 51]. Further, as shown by Stieltjes [50, 51], the solutions x\,x 2 , ■ ■ ■ ,xn 
of equations (3.2) are the roots of the N\h Hermite polynomial Hn(x). This highlights the connection between point 
patterns corresponding to a relative equilibrium and functions from classical mathematical physics. 
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3.2 Vortices of the same strength and mixed signs 

Now we consider the situation when the vortices at z\, z<i, ■ ■ ■ , z m have strength T > 0, i.e. there are m positive 
vortices, and the vortices at z m+ i, z m+2l ■ ■ ■ , z m+n have strength — T, i.e. there are n negative vortices, where m+n = 
N, so that 

r fe = ( r < for * = 1.2.->™ (3 .3) 
I— r, for k = m + 1, m + 2, . . . , m + n. 

We shall consider two important cases of equation (3.1) with these circulations. 

1 . Stationary vortex patterns that arise when dz*/di = 0, which we discuss in §3.2. 1 . 

2. Translating vortex patterns that arise when dz* /dt = v*, with v* a (complex) constant, which we discuss in 
§3.2.2. 

To study these vortex patterns we define the polynomials 

m n 

P(z) = H(z- Zj ), Q(z) = l[(z-( k ), (3.4) 

j=l k=l 

with (k = Zk+ m for k = 1, 2, . . . , n, which respectively have roots at z\, Z2, ... , z m , i.e. the locations of the positive 
vortices, and Q\, £2, . . . , Cn> i- e - the locations of the negative vortices. 

3.2.1 Stationary vortex patterns 

In this case, setting dz* /dt = in (3.1) gives 

X] — Z — = °' j = 1,2,..., m + n, 

fe=i z i Zk 

with Tfe given by (3.3). If P{z) and Q(z) are given by (3.4) then it can be shown that they satisfy the bilinear equation 

*^_ a ^ +J ,*O_0 «3,) 
dz 2 dz dz dz 2 

(for details see [8, 10, 32]). As noted in §2.2, polynomial solutions of equation (3.5) were studied by Burchnall 
and Chaundy [13], and later by Adler and Moser [3] — recall equation (2.8), see also Theorem 3.1. In the vortex 
dynamics literature equation (3.5) is sometimes known as Tkachenko's equation (cf. [8, 49]). Substituting P{z) = 
zi + cij-iz^ 1 + . . . + do and Q(z) = z k + bk-iz k ~ 1 + . . . + bo into equation (3.5) yields the leading order term 

- 1) - 2jk + k(k - l)]zi+ k - 2 = [(j - kf - (j + k)]z=+ k - 2 , 

thus j — k = n and j + k = n 2 and so j = ^n(n + 1) and k = ^n(n — 1), with n e Z. Thus the degrees of the 
polynomials P(z) and Q(z) are two successive triangular numbers. 

Burchnall and Chaundy [13] and Adler and Moser [3] showed that polynomial solutions of equation (3.5) are 
expressed in terms of the Adler-Moser polynomials, which were defined in §2.2, as given in the following result; see 
also [10, 12]. 

Theorem 3.1. Polynomial solutions of equation (3.5) are given by P(z) = 9 n +i(z; k„+i) and Q(z) — Q n (z; K n ), 
where 0„(z; K n ) is the Adler-Moser polynomial. 

Adler and Moser [3], using ideas on Darboux transformations due to Crum [23], also proved the following result, 
which led to the Wronskian representation of the Adler-Moser polynomials given in Definition 2.2. 

Theorem 3.2. The Schrddinger equation 

-^-4+u(z)^ = 0, (3.6) 
dz 

d 2 

with potential u(z) = —2 — j m 6n( z ; K n), where 8„(z; k„) is the Adler-Moser polynomial, has the eigenstate 
dz 

il>(z) = <d n+1 (z;K n+1 )/<d n (z;K n ). (3.7) 

We remark that substituting 

u(z) = -2^1nQ(z), #*) = §Jy. 
into the Schrddinger equation (3.6) yields the bilinear equation (3.5). 



7 



3.2.2 Translating vortex patterns 

In this case, setting dz* /dt — v*, with v* a constant, in (3.1) yields 

m+n 

i w r fc , 

— > =v , j = 1,2,..., m + n, 

27T1 f-J Zj - Zfe 

fe=l J 

with given by (3.3). If P(z) and Q(z) are given by (3.4) then it can be shown that they satisfy the bilinear equation 



dz* Q - 2 ^ + P J? + ^ \-z Q - P ^z)=^ (3 ' 8) 

where fj, = 2iriv* /V (for details see [8, 10, 32]). Substituting P{z) = z j + a j _ 1 z : >- 1 + . . . + a and Q(z) = 
z k + bk-iz k ~ 1 + . . . + b into equation (3.8) yields the leading order term 2 / u(j — fc)z- J+fe_1 and so j = k, i.e. 
the polynomials P(z) and Q(z) have the same degree. Theorem 3.3 shows that the degrees of the polynomials P(z) 
and Q(z) are triangular numbers, though it appears not to be understood why this restriction holds, independent of 
knowing the solution, in contrast to the situation for stationary vortex patterns in §3.2.1. 

Adler and Moser [3] proved the following result concerning polynomial solutions of equation (3.8); see also 
[10, 12, 39]. 

Theorem 3.3. Polynomial solutions of equation (3.8) are given by P(z) = Q n (z; fx, K n ) and Q(z) — Q n (z;n n ), 
where 6 rl (z; K n ) is the Adler-Moser polynomial and & n (z; \i, K n ) is the modified Adler-Moser polynomial. Further 
Schrddinger equation 

1 2 , , 

- -rin . | (■ = —/■- ■ 

dz 



-^4+u(^ = -M 2 V, (3-9) 



with potential u(z) = —2 — j m ®n{z; K n ), where Q n (z; n n ) is the Adler-Moser polynomial, has the eigenstate 
dz 

ip(z) = Q n {z; li, K n )/@ n (z; k„), (3.10) 

where & n (z; fi, K n ) is the modified Adler-Moser polynomial. 
We remark that substituting 

u { z) = -2-^ Q{ z), ^) = -^ r , 
into the Schrodinger equation (3.9) yields equation (3.8). 

3.3 Quadrupole Background Flow 

The equations of motion for N point vortices with circulations Tj at positions Zj, j = 1, 2, . . . , N in a background 
flow w(z) are 



dt 27ri^- Zj -z k 27ri ' 3 i,A---,iv. 

k=l J 



(3.11) 



Now suppose that dz* /dt = 0, i.e. stationary flow, w(z) — v*z*, with v* a (complex) constant, i.e. quadrupole 
background flow, and f^ is given by (3.3). In this case, setting v = LiT gives 

m+n 

E, Tk 
= -Tp.Zj, j = 1,2,..., m + n, (3.12) 
9!j — 

k=l 

and so using (3.3) yields 



m n 

E< 1 v-^ 1 
=-fiZj+y -, j = l,2,...,m, 
7- — ?t 7- — In 



(3.13a) 



-±-=^+f^—L-, j = l,2,...,n. (3.13b) 

vj t ._ 1 y Zfc 



<=i ,J ' fc=i 
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If P(z) and Q(z) are given by (3.4) then 



nz ^— ' y. — 7. a (\7. ^— ' z — 



dz 

d 2 P 

dz 2 

d 2 Q 




n n n 1 / m 1 



using (3.13) and so 



dz dz | ^^-^ ^ * - Cfc 



2 - Gc 



= 4^ - ^ ( P ¥) + - n)PQ, 
dz dz \dz dz J 

and so the polynomials P(z) and Q(z) satisfy the bilinear equation 

d2P Q - 2^ + P^ + 2,z f^Q - pW) - M m - n)PQ = 0. (3.14) 



dz 2 dz dz dz 2 \ dz dz 

Substituting P(z) = z rn + a m - 1 z m ~ 1 + . .. + a and Q(z) = z n + b n _ 1 z n ~ 1 + . . . + b into equation (3. 14) shows that 
the leading order term is identically satisfied and so the degrees of the polynomials P(z) and Q(z) are not restricted. 

Kadtke and Campbell [32] obtained some polynomial solutions of equation (3.14) in the case when m = n = 2 
and m = n = 4, though they also claimed that there were no solutions when m = n = 6. However, we have found 6 
pairs of solutions of equation (3.14) in the case when m = n = 6 and 12 pairs of solutions when m = n = 8 (for this 
case 2 pairs of solutions are given in [32]). The roots of these polynomials are invariant under reflections in the real 
and imaginary z-axes when m = n. In Tables 3.1 and 3.2 some polynomial solutions of equation (3.14) for fi = — | 
in the cases when m = n and m = n + 1 are given, respectively. All of these polynomials are Wronskians of the 
scaled Hermite polynomial, in fact they have the form (3.23), with (z) replaced by (z/y/2) and scaled so that 
the polynomials are monic. 

We note that equation (3.14) for \i = — 1 (which we can assume without loss of generality, by rescaling z if 
necessary), i.e. 

d2P Q - 2^ + P^-2z f^Q - pW) + 2(m - n)PQ = 0, (3.15) 



dz 2 dz dz dz 2 \ dz dz 

can be written in the form 

[D 2 - 2zD z +2(m-n)] P»Q = 0, (3.16) 

with D z the Hirota operator (2.12). Hence from (2.11) and (2.14), we can obtain three sets of solutions of equation 
(3.15), or equivalently (3.16), in terms of the generalized Hermite polynomials Hj t k(z) and the generalized Okamoto 
polynomials Q,j^{z), respectively defined in Definitions 2.3 and 2.4, which are given in Table 3.3. As shown below, 
it turns out that only two of these sets have the form (3.23). Further, Filipuk and Clarkson [26] show that there are 
solutions of equation (3.15), or equivalently (3.16), in terms of the symmetric Hermite polynomials Hk n (z) and the 
symmetric Okamoto polynomials Qk n (z), which are defined in Definitions 2.5 and 2.6, respectively. 

3.3.1 Relationship to the Schrodinger equation 

In §3.2 it was shown that stationary and translating vortex equilibrium solutions can be expressed in terms of the roots 
of the Adler-Moser and modified Adler-Moser polynomials and these were related to eigenstates of the Schrodinger 
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P(z) 


Q{z) 


m 


= n = 


2 


T 1 

z 2 + 1 


z 1 -! 


m 


= n = 


3 


z 3 + 3z 


z 3 










z 3 


z 3 - 3z 




m 


= n = 


4 


z 4 + 6z 2 + 3 


z 4 + 2z 2 - 1 










z 4 + 2z 2 - 1 


z 4 - 2z 2 - 1 










z 4 - 2z 2 - 1 


z 4 - 6z 2 + 3 




m 


= n = 


5 


z 5 + 10z 3 + 15z 


z 5 + 5z 3 










z 5 + 5z 3 


z 5 - 5z 










z 5 - 5z 


z 5 - 5z 3 










z 5 - 5z 3 


z 5 - 10z 3 + 15z 










z 5 + 2z 3 + 3z 


z 5 - 2z 3 + 3z 




m 


= n = 


6 


z b + 15z 4 + 45z z + 15 


z b + 9z 4 + 9z 2 - 3 










z 6 + 9z 4 + 9z 2 - 3 


z 6 + 3z 4 - 9z 2 - 3 










z 6 + 3z 4 - 9z 2 - 3 


z 6 - 3z 4 - 9z 2 + 3 










z 6 - 3z 4 - 9z 2 + 3 


z 6 - 9z 4 + 9z 2 + 3 










z 6 - 9z 4 + 9z 2 + 3 


z 6 -15z 4 + 45z 2 -15 










z 6 + 3z 4 + 9z 2 - 9 


z 6 - 3z 4 + 9z 2 + 9 




m 


= n = 


7 


z' + 21z b + 105z 3 + 105z 


z> + 14z b + 35z 3 










z 7 + 14z 5 + 35z 3 


z 7 + 7z 5 - 7z 3 - 21z 










z 7 + 7z 5 - 7z 3 - 21z 


z 7 - 21z 3 










z 7 - 21z 3 


z 7 - 7z 5 - 7z 3 + 21z 










z 7 - 7z 5 - 7z 3 + 21z 


z 7 - 14z 5 + 35z 3 










z 7 -14z 5 + 35z 3 


z 7 - 21z 5 + 105z 3 - 105z 










z 7 + 9z 5 + 15z 3 + 15z 


z 7 + 3z 5 - 3z 3 + 3z 










z 7 + 6z 5 + 15z 3 


z 7 - z 5 + 5z 3 + 15z 










z 7 + z 5 + 5z 3 - 15z 


z 7 - 6z 5 + 15z 3 




m 


= n = 


8 


z 8 + 28z b + 210z 4 + 420z 2 + 105 


z 8 + 20z b + 90z 4 + 60z^ - 


15 








z 8 + 20z 6 + 90z 4 + 60z 2 - 15 


z 8 + 12z 6 + 10z 4 - 60z 2 - 


15 








z 8 + 12z 6 + 10z 4 - 60z 2 - 15 


z 8 + 4z 6 - 30z 4 - 36z 2 + 9 








z 8 + 4z 6 - 30z 4 - 36z 2 + 9 


z 8 - 4z 6 - 30z 4 + 36z 2 + 9 








z 8 - 4z 6 - 30z 4 + 36z 2 + 9 


z 8 - 12z 6 + 10z 4 + 60z 2 - 


15 








z 8 -12z 6 + 10z 4 + 60z 2 -15 


z 8 - 20z 6 + 90z 4 - 60z 2 - 


15 








z 8 - 20z 6 + 90z 4 - 60z 2 - 15 


z 8 - 28z 6 + 210z 4 - 420z 2 


+ 105 








z 8 + 10z 6 + 30z 4 + 30z 2 - 15 


z 8 + 2z 6 + 30z 2 + 15 










z 8 + 4z 6 + 6z 4 - 12z 2 - 3 


z 8 - 4z 6 + 6z 4 + 12z 2 - 3 










z 8 - 2z 6 - 30z 2 + 15 


z 8 - 10z 6 + 30z 4 - 30z 2 - 


15 








z 8 + 8z 6 + 30z 4 + 45 


z 8 + 10z 4 - 15 










z 8 + 10z 4 - 15 


z 8 - 8z 6 + 30z 4 + 45 





Table 3.1: Some polynomial solutions of equation (3.14) for ji = — | and m = n 
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P(z) 


Q(z) 


m = 3, n = 2 


z 3 + 3z 


z z + 1 


m = 4, n = 3 


z 4 + 3 

z 4 + 6z 2 + 3 


z 3 - 3z 
z 3 + 3z 


m = 5, n = 4 


z 5 + 10z 3 + 15z 
z 5 - 2z 3 + 3z 
z 5 + 2z 3 + 3z 


z 4 + 6z 2 + 3 
z 4 - 6z 2 + 3 
z 4 - 2z 2 - 1 


m = 6, n = 5 


z° + 3z 4 + 9z 2 - 9 
z 6 + 5z 4 + 5z 2 + 5 
z 6 + 15z 4 + 45z 2 + 15 
z 6 - 5z 4 + 5z 2 - 5 


z° - 2z 3 + 3z 

z 5 - 5z 

z 5 + 10z 3 + 15z 
z 5 - 10z 3 + 15z 

z 5 - 5z 3 


m = 7, n = 6 


z 7 + 3z 5 - 3z 3 + 3z 
z 7 + 9z 5 + 15z 3 + 15z 
z 7 - 3z 5 - 3z 3 - 3z 
z 7 - 9z 5 + 15z 3 - 15z 
z 7 + 21z 5 + 105z 3 + 105z 
z 7 + z 5 + 5z 3 - 152 


z 6 - 3z 4 - 9z 2 + 3 
z 6 + 3z 4 - 9z 2 - 3 
z 6 - 9z 4 + 9z 2 + 3 
z 6 - 15z 4 + 45z 2 - 15 
z 6 + 15z 4 + 45z 2 + 15 
z 6 - 5z 4 + 5z 2 - 5 


m = 8, n = 7 


z 8 + 28z e + 210z 4 + 420z 2 + 105 
z 8 + 14z 6 + 42z 4 + 42z 2 + 21 
z 8 + 4z 6 + 6z 4 - 12z 2 - 3 
z 8 + 10z 6 + 30z 4 + 30z 2 - 15 
z 8 - 14z 6 + 42z 4 - 42z 2 + 21 
z 8 + 8z 6 + 30z 4 + 45 
z 8 - 2z 6 - 30z 2 + 15 
z 8 - 14z 4 - 7 
z 8 + 7z 6 
z 8 - 7z 6 


z 7 + 21z 5 + 105z 3 + 105z 
z 7 + 7z 5 - 7z 3 - 21z 
z 7 - 3z 5 - 3z 3 - 3z 
z 7 + 3z 5 - 3z 3 + 3z 
z 7 - 21z 5 + 105z 3 - 105z 
z 7 + z 5 + 5z 3 - 15z 
z 7 - 9z 5 + 15z 3 - 15z 
z 7 - 7z 5 - 7z 3 + 21z 
z 7 - 21z 3 
z 7 - 14z 5 + 35z 3 



Table 3.2: Some polynomial solutions of equation (3.14) for /i = — \ and m = n + 1 



P(z) 


Q(z) 


m — n 


H j,k( z ) 
H jik +i(z) 
Hj. k +i(z) 


Hj+i,ki z ) 

Hj,k{z) 
H j+lt k{z) 


-k 

i 

j-k 


%+i,fc( 2 ) 
&j,k(z) 


^j.k(z) 
%,fe(z) 
Qj+i,k+i(z) 


j -2k + 2 
-2j + k + l 
-J -k — 2 



Table 3.3: Solutions of (3.15) and (3.16) in terms of generalized Hermite and generalized Okamoto polynomials. 
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equations (3.6) and (3.9). Now in order to describe quadrupole background flow equilibrium solutions we shall relate 
solutions of equation (3.15) to eigenstates of a Schrodinger equation. 
Consider the Schrodinger operator 

C = -^+u(z), (3.17) 
dz 

with a potential u(z) which is meromorphic in C. The Schrodinger operator £ (3.17) is said to have trivial monodromy 
if all the solutions of the corresponding Schrodinger equation 

d? 



Ctp = — + u{z)tp = Xip, (3.18) 



are also meromorphic in C for all A. Further, such an operator is said to be monodromy-free. 

Duistermaat and Griinbaum [25] considered the problem of the classification of Schrodinger operators for rational 
potentials decaying at infinity, i.e. (3.17) with limi^^ u(z) = 0, and obtained the following result. 

Theorem 3.4. Every Schrodinger operator C with trivial monodromy for rational potentials decaying at infinity has 
the form 

^-^ + E^ (3-19) 

dz 2 f-f {z-zjf 

where 1/1,1/2, ... ,Vn are integers. Further all such potentials are obtained by finitely many Darboux transformations 
applied to the zero potential. 

Oblomkov [44] generalized Theorem 3.4 to the case of rational potentials with quadratic growth at infinity and 
obtained the following result. 

Theorem 3.5. Every Schrodinger operator C with trivial monodromy, and with a quadratically increasing rational 
potential, i.e. u{z) = z 2 + R(z), with lim^i^^ R(z) = 0, has the form 

C = --^2 +z 2 -2^ lnW(H kl (z),H k2 (z), H k( (z)) , (3.20) 
dz dz 

where k\ , ki, ■ ■ ■ , ke are a sequence of distinct positive integers and H k (z) is the kth Hermite polynomial. Further all 
such potentials are obtained obtained by finitely many Darboux transformations applied to the potential u = z 2 . 

Then, using Theorem 3.5 and a corollary of results on Darboux transformations due to Crum [23] we have the 
following result (see also Veselov [54]). 



Theorem 3.6. The Schrodinger equation (3.18) with potential 

dz 2 

where k\, ki, ■ ■ ■ , ke are a sequence of distinct positive integers and H k (z) is the kth Hermite polynomial, has the 
eigenstates 



u(z) = z 2 - 2— In W (H kl (z),H k2 (z),..., H ke (z)) , (3.21) 



{z) W(H kl (z),H k2 (z),...,H ki (z),H ke+1 (z)) 
nZ) W(H kl (z),H k2 (z),...,H ke (z)) ° XP[ 2 ">' (3 ' 22a) 



^ W(H kA z),H k2{ z),...,H ke{ z)) 6XP (** ) ' (3 ' 22b) 

with kg + i a positive integer such that k^ + i ^ kj, for j — 1,2, ... ,£, for the eigenvalues A = 1 + 2(ki +1 — £) and 
A = 2(£ — k e _i) — 1, respectively. 

It should be noted that the potential (3.21) may possess singularities on the real axis, so that the terms "eigenstate"' 
and "eigenvalue" should not be understood in quantum mechanics sense. The associated regularity conditions have 
been studied by Adler [5] and Dubov, Eleonskii, and Kulagin [24]. 

We remark that substituting 

u(z) = z 2 2^- 2 lnQ(z), = cxpH, 2 ), 

into the Schrodinger equation (3.18) yields equation (3.15) with m — n = A — 1. Consequently we have the following 
result (see also Loutsenko [38, Proposition 4]). 



12 



Theorem 3.7. The bilinear equation (3.15) has polynomial solutions in the form 

P(z) = W (H kl (z), H k2 (z), . . . , H kt (z),H ke+1 (z)) , (3.23a) 
Q(z) = W (H kl (z), H k2 {z), ...,H ke (z)) , (3.23b) 

where ki, k 2 , . . . , kg, kg + i are a sequence of distinct positive integers and H k (z) is the kth Hermite polynomial. 
Further the degrees of the polynomials P(z) and Q(z), respectively m and n, are given by 

i+i i 
m = j>-^(*+l), n = J>-§*(*-l) 

3 = 1 3=1 

i.e. in — n = ke+i — t 

We note that the polynomials P(z) and Q(z) given by (3.23) do not involve arbitrary parameters, in contrast to the 
Adler-Moser polynomials which give polynomial solutions of equations (3.5) and (3.8). 

Theorem 3.7 would appear to give a complete classification of polynomial solutions of the bilinear equation 
(3.15). As remarked above, all the polynomials in Tables 3.1 and 3.2 have the form (3.23), with H kj (z) replaced 
by H kj (z/a/2) and scaled so that the polynomials are monic. 

However there are additional polynomial solutions of equation (3.15) in terms of the generalized Hermite polyno- 
mials H klik2 (z) and the generalized Okamoto polynomials Q kl >k2 (z). In Table 3.3, there are three pairs of generalized 
Hermite polynomials and three pairs of generalized Okamoto polynomials which satisfy equation (3.15), though only 
two pairs of each are of the form (3.23) in which the wronskian defining P(z) has one more Hermite polynomial than 
the wronskian defining Q(z). Specifically the polynomials 

P{z) = H klM (z) = W (H kl , H kl+1 , tf fcl+fca _i) , (3.24a) 
Q(z) = H kl+1M (z) = W (H kl+1 ,H kl+2 , H kl+k2 ) , (3.24b) 

where the wronskians defining P(z) and Q(z) have the same number of Hermite polynomials, and 

P{z) = fl klM (z) = W{H U H 4 , H 3kl _ 2 ,H 2 , H 5 , flafa-i), (3.25a) 
Q(z) = n kl+ i ik2+ i(z) = W(Hi, H4,, . . . , i?3fe 1+ i, H 2 , H 5} . . . , H 3k2+2 ), (3.25b) 

where the wronskian defining P(z) has two fewer Hermite polynomials than the wronskian defining Q(z), are pairs 
of solutions of equation (3.15), yet are clearly not of the form (3.23). An analogous result also holds for the symmetric 
Hermite polynomials H/~ 2n (z) and the symmetric Okamoto polynomials Q kn (z). Specifically the polynomials 

P(z) = H klMt ... Mn {z) 

= w[{H kl+j (z)}^-\{H k 

i+k 2 + ... + k 2n -i+j 

(z)}^ Q , (3.26a) 

Q(z) = H kl+1MM+ i M ,„ Mn _ 1+ i Mn (z) 

= W ({H kl+j {z)}% v {H kl+k2+ka+3 (z)} k / =v {H kl+k2+ ... +k2n _ 1+j (z)}*l\) , (3.26b) 

where the wronskians defining P(z) and Q(z) have the same number of Hermite polynomials, and 

P(z) = £lk 1 ,k2,...,k n (z) 

(n+l)j+n 

Wl-Io 1 )' (3 ' 27a) 

Q(z) = 0fe 1 +i,fe 2 +i,...,fc„+i(2) 

= W ({-ff(n+l)j+l(^)}^l 0) {#(n+l)j+2(z)}*l , • • • , { H {n+\)j+n{z)} k ^ , (3.27b) 

where the wronskian defining P(z) has n fewer Hermite polynomials than the wronskian defining Q(z), are pairs of 
solutions of equation (3.15), yet also are not of the form (3.23). Hence we have demonstrated that there are additional 
polynomial solutions of equation (3.15) to those given in Theorem 3.7. The natural question is whether there are any 
more? 

Using MAPLE, we have also found the polynomial solutions of equation (3.15) with m = n and with m = n + 1, 
given in Tables 3.4 and 3.5, respectively, which have the property that either the polynomials P(z) or Q(z), or both, 
have multiple roots, other than z = 0, and the polynomials P(z) and Q(z) have common roots, other than z = 0. 
Further these polynomials do not appear to be wronskians of Hermite polynomials. 
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P(z) 


Q(z) 


(z 2 + V5z + 2) 3 
(z 2 + | + iiVT5) 3 

(z + ±V6) 6 
(z 2 + l) 3 (z 2 -l) 
z 3 (z 2 + |) 3 

z(z 2 + | + iiV3T) 3 
(z 2 + l) 3 (z 2 -l) 

z 3 [z 2 + |(l + iV7)] 3 
(z 2 + |)6 

(Z 2 + | )6(Z 2_| )3 


(z 2 + \/5z + 2)(z 4 + 2\/5z 3 + 6z 2 - 3) 
(z 2 + | + ii^)[z 4 -(|-ii^)z 2 + |] 
(z 3 + |V6z 2 + |z-|V6)(z + iV6) 3 
(z 2 + l)(z 2 -l) 3 

z(z 2 + | )(z 6_| z 4 + | z 2_| ) 

z 3 (z 2 + | + iiV3T)(z 2 -l + iiV3T) 
(z 2 + l)(z 2 -l) 3 
z(z 2 - §) [z 2 + 1(1 + iV7)] [z 4 - |(1 - iV7)z 2 - |] 

(z2+ 3 ) 3 (z 6_3 z 4 + fz2 _| ) 
(z 2 + | )3(z2 _| )6 



Table 3.4: Some additional polynomial solutions of equation (3.15) with m = n. 



P(z) 


Q{z) 


(z + 1) 3 

(z 3 + 2V3z 2 + fz + V3)(z + V3) 
(z 2 + I) 3 
(z 2 + VT0z + 3) 3 

(z + \VUf 

z(z 2 -f)(z 4 + f) 
(z + 2) 6 (z + 3) 
(z-2) 6 (z-3) 


(z + 2)(z + l) 
(z + V3) 3 
z(z 2 + I)(z 2 -§) 
(z 2 + VWz + 3)(z 2 + fVTOz + 6)(z + i^/To) 
(z 2 + f\/l4z + 8)(z + ±\/l4) 3 
(z 2 -§) 3 
(z + 2) 3 (z + 3) 3 
(z-2) 3 (z-3) 3 



Table 3.5: Some additional polynomial solutions of equation (3.15) with m = n + 1. 
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4 Discussion 



In paper we have investigated the relationship between equilibria of vortex patterns and the roots of polynomials 
associated with rational solutions of the Painleve and soliton equations. The examples which we have considered in 
this paper illustrate the power of the "method of polynomials" when it can be used. The polynomials with roots at 
the positions of the vortices, are obtained from the basic vortex equilibrium equations and then their roots give the 
locations of the vortices. 

When all the vortex circulations have the same absolute value, stationary vortex configurations and uniformly 
translating vortex configurations are known to be given by the zeros of the Adler-Moser polynomials and related to 
eigenstates of the Schrodinger equation for rational potentials decaying at infinity. 

For quadrupole background flow and all the vortex circulations have the same absolute value, vortex configurations 
were shown to be given by the zeros of polynomials expressed as wronskians of Hermite polynomials and related to 
eigenstates of the Schrodinger equation for rational potentials with quadratic growth at infinity. Further using results on 
the generalized Hermite polynomials H m . n (z) and generalized Okamoto polynomials fl m ,n(z), which are associated 
with rational solutions of Piv (1-3), it was shown that there are additional solutions than might be anticipated from 
earlier work of Loutsenko [38], Oblomkov [44], and Veselov [54]. As mentioned above, a natural question is whether 
there are any more? 

An interesting open question is whether there is a relationship between vortex dynamics and polynomials associ- 
ated with rational solutions of other soliton equations, such as the Boussinesq equation (1.4) and the NLS equation 
(1.5). Since P IV (1-3) arises as a similarity reduction of the Boussinesq and NLS equations, then one can express ratio- 
nal solutions of these equations in terms of the generalized Hermite and generalized Okamoto polynomials. However 
both the Boussinesq and NLS equations also have rational solutions which involve arbitrary parameters, analogous to 
the Adler-Moser polynomials which describe rational solutions of the KdV equation (1.1). For the NLS equation (1.5) 
an extension of the generalized Hermite polynomials, see [20], and for Boussinesq equation (1.4) these polynomials 
involve an extension of the generalized Okamoto polynomials, see [21]. These are obtained by replacing the Hermite 
polynomial H m (z), which is defined by the generating function 

g HrrM>r = _ 

^— mi 

m=0 

by the polynomial $ m (z) defined by 

m=0 ' y j=3 

where Kj, for j > 3, are arbitrary constants in the definitions of the generalized Hermite and generalized Okamoto 
polynomials. (This is similar to the relationship between the Yablonskii-Vorob'ev polynomials and the Adler-Moser 
polynomials where the polynomial tp n (z) defined in equation (2.1) is replaced by the polynomial 9 n (z) defined in 
equation (2.5).) Therefore the structure of rational solutions of the Boussinesq equation (1.4) and the NLS equation 
(1.5) is different to those of the KdV equation (1.1) and so other vortex dynamics can possibly be described in terms 
of the associated generalized polynomials. 

Another interesting open question is whether there exist polynomial solutions of the equations 

fZ Q _ 2A 1£10 + A ^f0 + „ (Vq _ + KPQ . o, (4.3) 

dz dz dz dz \dz dz J 

where A ^ and n are constants, which respectively reduce to equations (3.5), (3.8), and (3.14), when A = 1. 
Loutsenko [39] has obtained some polynomial solutions of equations (4.1) and (4.2) in the case when A = 2, which 
corresponds to the case of vortices with two strengths of unequal magnitudes. 

Patterns of points corresponding to the equilibria of point vortices reveal interesting connections to the classical 
theory of polynomials and to some of the special polynomials that arise in soliton theory. While some of these 
connections have been discovered, we believe there is much more that awaits to be found given the large number of 
results which are known in soliton theory. 
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